Inverse energy cascade in three-dimensional isotropic turbulence 
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In turbulent flows kinetic energy is spread by nonlinear interactions over a broad range of scales. 
Energy transfer may proceed either toward small scales or in the reverse direction. The latter case is 
peculiar of two-dimensional (2D) flows. Interestingly, a reversal of the energy flux is observed also in 
three-dimensional (3D) geophysical flows under rotation and/or confined in thin layers. The question 
is whether this phenomenon is enforced solely by external anisotropic mechanisms or it is intimately 
embedded in the Navier- Stokes (NS) equations. Here we show that an inverse energy cascade occurs 
also in 3D isotropic flow. The flow is obtained from a suitable surgery of the NS equations, keeping 
only triadic interactions between sign-defined helical modes, preserving homogeneity and isotropy 
and breaking reflection invariance. Our findings highlight the role played by helicity in the energy 
transfer process and show that both 2D and 3D properties naturally coexist in all flows in nature. 



Inviscid invariants of the NS equations are crucial in 
determining the direction of the turbulent energy transfer 
PQ. In some cases, as for fully isotropic and homogeneous 
turbulence in 2D the presence of two positive-defined in- 
variants (energy and enstrophy) does not allow a sta- 
tionary transfer of both quantities, neither to small nor 
to large scales [2 j. In presence of two fluxes, they must 
necessarily flow in opposite directions [3H7] and this re- 
mains true even for turbulent systems in non-integer di- 
mensions obtained by fractal Fourier decimation [8 . The 
fluid equations possess two inviscid invariants also in 3D: 
energy and helicity (i.e. the scalar product of velocity 
and vorticity). The inviscid conservation of helicity was 
discovered relatively recently p~6j [17] . At variance with 
energy, helicity is not positive denned. This allows for a 
simultaneous small-scale transfer of energy and helicity, 
as confirmed by results of two-point closures [17-19 and 
direct numerical simulations [20 , 21]. Nevertheless, a re- 
versal of the flux of energy has been observed in geophysi- 
cal flows subject to earth rotation [9j[10] as well as in shal- 
low fluid layers [TTHT5] . In both cases, this phenomenon 
is accompanied by strong anisotropic effects and by a 
substantial two-dimensionalization of the flow, induced 
either by the rotation or by the effects of confinement. 
Moreover, rotations injects fluctuations in the helical sec- 
tor while a perfect two-dimensional flow has vanishing 
point-wise helicity, being vorticity always ortogonal to 
velocity. The role played by helicity in the energy trans- 
fer mechanism of 3D flows has attracted a broad scientific 
interest (see, e.g., [21 and reference therein). Dynamical 
systems have been developed to study in details energy 
and helicity transfer at high Reynolds numbers [22j [23] . 
Further, speculation connecting the existence of intermit- 
tent burst in the energy cascade induced by "local" he- 
licity blocking mechanism have been proposed [22 . De- 
spite these important contributions, the understanding 
of the phenomenology of helicity remains "mysterious" , 



as summarized in the conclusion of a recent state-of-the- 
art numerical study [21 . Here we present theoretical 
and numerical evidences of a new phenomenon induced 
by helicity conservation: a statistically stationary back- 
ward energy transfer can be sustained even in 3D fully 
isotropic turbulence. 

The starting point of our analysis is the well-known 
helical decomposition [19 of the velocity field v(x), ex- 
panded in Fourier series, u(k): 

u(k) = u+(k)h+(k) + u-(k)h~(k) (1) 

where /i ± are the eigenvectors of the curl operator ik x 
/i ± = ^zkh^. In particular, we choose /i ± = x k ± iO, 
where is an arbitrary versor orthogonal to k which sat- 
isfies the relation v(k) = —u(-k) (necessary to ensure 
the reality of the velocity field) . Such requirement is sat- 
isfied e.g. by the choice v = z x k/\\z x k\\ : with z an 
arbitrary vector. In terms of this exact decomposition of 
each Fourier mode energy, E = J d 3 x\v(x)\ 2 , and helic- 
ity, H = J d 3 xv • w, where w is the vorticity, are written 
as: 

E = J2 k \u + (k)\ 2 + \u-(k)\i; 
H = Z k k(\u+(k)\*-\u-(k)n 1 ' 

Similarly, the non-linear term of the NS equations can be 
exactly decomposed in 4 independent classes of triadic in- 
teractions, determined by the helical content of the com- 
plex amplitudes, u Sk (k) with = =b (see [19 ]). Among 
three generic interacting modes u Sk (k),u Sp (p),u Sq (q), 
one can identify 8 different helical combinations (s& = 
±, s p = =b, s q = zb). Among them, only four are indepen- 
dent because of the symmetry that allows to change all 
signs of helicity simultaneously. Let us now consider the 
dynamics of an incompressible flow V • v = 0, which is 
determined by a decimated NS equation in which all in- 
teractions between modes have been switched off except 
for those with a well defined sign of helicity, e.g. positive 



2 




VORTICITY HELICITY 




FORWARD ENERGY TRANSFER 



FIG. 1: Comparison between helicity and vorticity fields for normal NS turbulence (bottom row) and for inverse cascade 3D 
turbulence (top row). We also represent a pictorial scheme of triads responsible for the inverse cascade regime as suggested 
in [19]. Inverse 3D cascade is possible with triads where all helical components have the same sign. In this case the middle 
wave number (here u + (p)) transfers energy also to smaller wavenumbers (u + (k)), at difference from what happens for triads 
with different helical components where it is always the smallest wavenumber u + (k) that transfer energy to higher modes 
u~(q),u + (p). 
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FIG. 2: (color online) Non-stationary spectrum in the in- 
verse energy cascade regime. Red, dashed line represent k~ 5 ^ 3 
slope. 



+). We define the projector on 



positive/negative helicity states as 



(3) 



where • stands for complex conjugate. Then we project 
the velocity field into its positive helicity component: 



v + (x) = ^e ikx V + u(k)', 



and we consider the decimated NS equations: 

d t v+ = (-v+ • Vv + - Vp) + + vAv+ + /+ 



(4) 



(5) 



where v is the viscosity, p is the pressure and / is the 
external forcing stirring the fluid around a wavevector 
kf. The non-linear term and the forcing are projected on 
the positive helicity states with the same procedure fol- 
lowed for the velocity field Q. The resulting system has 
two positive definite invariants, see Eq. the energy 
and the helicity, H = J2 k k\u+(k)\ 2 , and contains only 
interactions between positive helicity modes. Helicity be- 
comes a coercitive quantity: the decimated NS equations 
cannot sustain a simultaneous forward cascade of energy 
and helicity, for the same arguments which forbid the 
existence of a simultaneous forward cascade of energy 
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and enstrophy in 2D turbulence [2j[19]- Therefore, the 
dynamics of Eq. ([5| should display a double cascade phe- 
nomenology, characterized by an inverse energy cascade 
with Kolmogorov spectrum E(k) ~ /c -5 / 3 for k <C fc/, 
and a direct helicity cascade with a k~ 7 ^ 3 spectrum for 
k ^> kf. It is interesting to note that, at variance with 
usual 3D NS dynamics, such flow should not display dissi- 
pative anomaly for kinetic energy, i.e. energy dissipation 
should vanish in the limit v — » 0. Indeed, the direct he- 
licity cascade carries also a residual, non-constant flux of 
kinetic energy toward small scales which decays as k~ x 
and therefore vanishes in the high Reynolds number limit. 
As a consequence, one may speculate that the decimated 
NS equations posses a less singular spatio-temporal evo- 
lution. Numerical simulations has been performed with 
a fully-dealiased, pseudo-spectral code at resolution 512 3 
on a triply periodic cubic domain of size L = 2tt. The 
flow is sustained by a random Gaussian forcing, with 
(fi(k,t)fj(q,t')) = F(k)S(k - q)S(t - t)Q itj (k), where 
Qij(k) is a projector assuring incompressibility and F(k) 
has support only in the high wavenumber range \k\ G 

[kmin = 25 '. kfnax = 32]. 

A visual inspection of the helicity and vorticity fields 
offered in Fig. [I] shows the differences between the for- 
ward cascade, which develops in standard 3D NS equa- 
tions forced at large scale, and the novel 3D inverse 
cascade regime obtained from the decimated NS equa- 
tion ([5| forced at small scales (see Fig. [I]). The latter 
does not posses any filamentary structure in the vortic- 
ity field, witnessing the fact that the vortex stretching 
mechanism, which is responsible for the forward cascade 
in standard 3D systems, is here reduced. Similarly, the 
forward regime does not possesses any coherent helicity 
signal at variance with the inverse regime which shows 
strong non-homogeneity in the helical spatial distribu- 
tion. 

In Fig. [2] we show a typical evolution of the energy 
spectrum obtained from Eq. ([5| by initializing the flow 
with energy only at high wave-numbers. The develop- 
ment of an inverse cascade with a Kolmogorov spectrum 
E(k) ~ k~ 5 / 3 is unambiguous. 

In absence of a large-scale dissipative mechanism, the 
inverse cascade would accumulate the kinetic energy 
in the lowest available mode, originating a condensed 
state [14] . In order to avoid this phenomenon we made a 
second series of numerical simulations, adding an hypo- 
viscosity at large scale ex A -1 v. In such a case, the total 
kinetic energy becomes stationary as shown in Fig. [3j and 
is equally distributed among the three velocity compo- 
nents, showing that the flow is fully isotropic. This allows 
to study scaling properties without having to cope with 
anisotropic sub-leading contributions [26 . In the inset 
of Fig. |3j we show the stationary energy flux in Fourier 
space, defined as II(fc) = (d/dt) E(p)dp where time 
derivative is computed by taking into account only the 
non-linear terms of Eq. ([5| . The negative plateaux in the 
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FIG. 3: (color online) Evolution of the three components of 
the turbulent kinetic energy as a function of time, {(vi) 2 ), with 
i = x (red, solid line) i = y (green, dashed line) and % — z 
(blue, dotted line). Inset: energy flux, II(fc), in the Fourier 
space. Notice the clear negative plateaux in the inert ial range 
k < k f . 



inertial range of wave-numbers is a clear indication of the 
large-scale energy transfer. 

The inverse cascade which arise from Eq. ([5| is not in- 
termittent. The probability distribution functions (pdfs) 
of the longitudinal velocity increments S r v = [v(x + r) — 
v(x)] • r at distance r within the inertial range are self- 
similar and almost Gaussian (see inset of Fig. [4|. The 
scaling of the second and the fourth order moment of 
velocity increments S^O") = ((^r^) 2 ); S^r) = ((S r v\^) 
follow the dimensional scaling S p (r) ~ r p / 3 (see Fig. kl. 
This is a signature of all known inverse cascades: when 
fluctuations are transferred from faster to slower degrees 
of freedom [27] . Previous studies have shown the possibil- 
ity to produce large scale motion by non-parity invariant 
small-scales forcing only at small Reynolds numbers or in 
the quasi-linear regime [25]. Conversely, our results do 
not trivially originate from the projection of the forcing 
on the positive helicity states but is a genuine effect of the 
non-linear dynamics. To assess this issue we performed 
a test simulation of the complete NS equation with the 
same projected forcing used in Eq. ([5|. After an initial 
transient in which part of the energy accumulates at the 
forcing scale, a direct cascade sets in and all the energy 
injected is transferred toward small scales. This excludes 
the possibility that the forcing alone could be responsible 
for the inverse energy transfer observed in the decimated 
NS equation. 

In conclusions, we have presented theoretical and nu- 
merical evidences that a screwed version of the NS equa- 
tions, such that only modes with a given sign of helicity 
are retained, displays inverse energy transfer mechanism. 
This phenomenon, which has been previously observed 
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only in 2D turbulence or in strongly anisotropic 3D flows 
under bidimensionalization effects, is here observed for 
the first time in a fully isotropic 3D system and is intrin- 
sically connected to the non-linear dynamics of all flows 
in nature. 

The scientific impact of our findings is manifold. First, 
it allows to highlight those backward events in the energy 
transfer mechanism which are known to exist also in un- 
truncated NS equations and that are one of the main the- 
oretical and applied challenges, see e.g. [24] for the case 
of sub-grid modelling in Large Eddy Simulations. Sec- 
ond, the link between backward energy events with the 
helical nature of triad interaction, shows the key role of 
the coupled energy-helicity dynamics. Third, by clearly 
detecting which triadic interaction is responsible for for- 
ward and backward energy transfer, we pave the road for 
closure and analytical approaches aimed at understand- 
ing the whole energy transfer distribution. 

This study also opens the way to further investigations. 
An obvious extension would be to integrate Eq. ([5| with 
a large scale forcing. In this case a pure forward helicity 
cascade must develop, provided that energy is removed at 
the forcing scale to avoid pile- up of fluctuations. More in- 
teresting, one could consider the case of a complementary 
decimation with respect to the one discussed here, i.e. 
eliminating only those triads that transfer energy back- 
ward. It is very tempting to speculate that such system 
could display a direct energy cascade with reduced -or 
even vanishing- intermittency, because one has removed 
all the obstacles, i.e. those events in which the forward 
energy transfer is stopped and/or reversed by the inter- 
action with the helicity flux. Numerical simulations ex- 
ploring these cases are ongoing and will be reported else- 
where. Finally, similar decomposition may shed lights 
also in the evolution of conducting fluids, where three 
invariants, kinetic plus magnetic energy, cross helicity 
and magnetic helicity are known to produce a reach phe- 
nomenology [28] . We acknowledge useful discussion with 
U. Frisch. L.B. acknowledge the kind hospitality from 
the Observatoire de la Cote d'Azur in Nice where part 
of this work was done. We acknowledge the European 
COST Action MP0806. 
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